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1 E-mail PESANDO@NBIVAX.NBI.DK, 22105::PESANDO, 31890::I_PESANDO. 
2 Work supported by the EU grant ERB4001GT930141. 



1 Introduction 



One of the most interesting and still not fully understood problem of Quan- 
tum Field Theory concerns the calculation of the mass of relativistic bound 
states. This problem is even more complicated to solve in the case of QCD, 
the SU (3) gauge theory of quarks and gluons, that describes the strong in- 
teraction physics where the basic constituents are confined and perturbation 
theory cannot be used. The expectation is that the low energy spectrum 
consists of colorless mesons and baryons, but up to now the only available 
method seems to be based on numerical calculation using the lattice theory 
formulation of the theory. 

The large- N expansion technique proposed by 't Hooft [1], several years 
ago seems to be the most promising approach to obtain analytic results on 
the hadron spectrum. It is well known that in the limit in which the num- 
ber N of colours becomes very large the theory becomes much simpler, in 
the sense that only planar Feynman graphs survive. Furthermore, when 
N — > oo the theory only contains colorless, stable and noninteracting mesons 
with two-body decay and scattering amplitudes proportional to h= and jj, 
respectively (there is also a way of describing baryons as solitons of the ef- 
fective Lagrangian in the large iV limit, but this will not interest us in this 
context). 

Unfortunately, all efforts to try to solve four-dimensional large- N QCD 
have failed up to now. The main problem is to find the semi-classical con- 
figuration, namely the master field ||, from which the action is dominated 
in the large- iV limit and whose fluctuations around the vacuum should give 
the particle spectrum of the theory. 

The main reason of this failure is due to the fact that no method has 
yet been found to solve matrix models for a space time dimension D > 
2. On the other side it is well known that the large N expansion can be 
explicitly performed in vector models as for instance the O(N) vector model 
(see for instance ref. J| and references therein), the two-dimensional CP^ -1 
model || and also in QCD 2 fH, || and |7], ||, |i| with matter in the fundamental 
representation of the gauge group Some results can be obtained with both 
fermionic and bosonic adjoint matter JlO| though the model is a real 2D 
matrix model. 

In particular this set of models are solved by means of two slightly differ- 
ent methods. The vector-like models are solved in the large N expansion by 
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introducing a local "colourless" composite field and by explicitly integrating 
over the fundamental fields (See for instance Refs.[|J and 0): with the help 
of the composite local field we can extract the explicit N dependence and 
apply the saddle point technique. 

The QCZ^-hke models are characterized by the fact that the number 
of components of the non abelian gauge field goes to infinity in the large 
N expansion and from this viewpoint they are matrix models. The gauge 
field in two dimensions has no physical degrees of freedom and therefore can 
be eliminated by using its classical equation of motion. In this way one 
gets a non local Coulomb interaction that is quartic in terms of the vector- 
like matter fields. The theory can then be solved in the large N limit by 
introducing a bilocal composite field as discussed in Refs. |7], |9|. We can 
rephrase this by saying that the kernel of the solubility lies in the fact that 
the matrices of the theory (the gauge fields) are expressible as tensor product 
of vectors. The same happens in the generalized QCD , even if the theory 
seems a 2 matrices model, it is actually a nonlocal vector model. 

In this paper we present a solution of the most general case of generalized 



QCD2 using the bilocal method (see |L3[ for an extensive application of the 
method to the other known cases). 
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2 The master field of the general QCD2:the 
simplest case. 

The usual QCD 2 action can be written in the first order formalism as [] 



S = [ d 2 x i—trJEF) - ^trJE 2 ) + ^ a (dj- m i lU a \ (2.1) 
J [ 8n An 1 j 

where F„ v = d^A v - d v A„ + i[A^A v \ } F = F^, if^ = ^(id„ t AB - 

A a T^ B ), a,b = 1...N 2 are the indices of the adjoint representation of the 
colour group, A,B = 1...N run over the fermionic representation of the 
colour U(N) and a, b... are flavour indices, and tr c is the color trace. 
From the previous form it can be directly generalized to ( ||1 1|| ) 

S = Jd 2 x ^tr c (EF)-^Y, fntr c (E n )+r(if-m l 11)^ j (2.2) 

To solve this theory we choose the usual gauge At = 0. and we can 
rewrite the previous expression as 



S=f d 2 x tr e (Ed-A+) - ^£ f n tr e (E») 

J {An An ^ 

+ ^) a (ii -m t )f - 4(f7_ff)l (2.3) 



1 Conventions. 



x ± = x T = -^{x° ± x 1 ) A^Bfj, = A B - AxBx = A + B_ + A + B_ 
a/2 



e 01 



V2\ ( \ /0 1\ / 1 

+ = ^ ) 7 = { V2 ) 70 = { 1 ) 71 = (, -1 

( 1 \ I ±15 

75 = -7071 = y _! ) Pr,l = —j- 



= / d 2 z 



V 



dp 
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where tr a and tr a are the traces over the flavour and spin indices. 
If we integrate over A+, we get the constraint 

d x - (— eA + ^ a 7 _r> a = (2.4) 

in the form of a delta function in the path integral. Since E is in the adjoint 
representation of U(N), the constraint ( |2.4| ) can be rewritten as 

d x - (J?- E AB ^j + i> Ba 1 _ip Aa = (2.5) 

Now we can integrate easily over E AB getting the effective action for the ip 
fields 

S eff = [d 2 xU Aa (ifi-m a )ij Aa 



- 77 Y,(%) n fn(dz^ A2ai i^ Aiai ) x ■ ■ • (d: V ia "7-^" a ") a 

(2.6) 

The interaction term suggests to introduce the (composite) bilocal field 

u m , v (x, y) = u PQ = « M ) = ( ) 

(2.7) 

where P = (xaa), Q = (yb/3), and 

PpQ(^) = P X aa,y bP (u) = (~^=) &{x ~ y)dZ\u, y) 6°* ( 7 -) a/3 (2.8) 



The effective action can then be rewritten as 
1 



^e// = j x ^\x-y)tr a , a ({ilf) x -m a )U{x,y)) + ^ 



(2.9) 



where Tr = tr P = tr x tr a tr a , and 

/ --^m a 5 ab S 2 (x-y) i5 ab d x -5\x-y) 
U-\\E>pq\\ - y i S ab d x+ 5 2 {x - y) -±m a 5 ah 5 2 {x-y) 



(2.10) 
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In order to get the effective action for the field U, we have to compute the 
jacobian of the transformation J[U] from ip to U. This can be accomplished 
as in Ref.s (0, \F3). The result is 



J[U) oc j [dM] exp {N[i Tr(MU) + Tr log M] } (2.11) 

The final form for the effective action is given as a functional of the U 
and M fields, and it is 

jjS eff =Tr(DU + ^-J2fn J[P{x)U] n + M U - i log Jtij (2.12) 

From this action we can immediately compute the saddle point equations: 
U M -i 31 = 

D + ^Zn nf n J x [P(x)ur- l P(x) +M = (2 ' 13) 



that yield 

2 

DU + ^-Y, n fnl (P(x)UT = -* B (2-14) 

47T ~ Jx 

Following the approach used in the usual case (f|), we introduce the self- 
energy £ = \ ( ^ L ~ ) through the definition 



R 



U = -i(D + Z)- 1 (2.15) 



o 



Inserting it in eq. ( |2.14| ) we get the following eq. for E 



n-l 



that yields immediately = T_ = = because of the 7_ included in 
P(x). Assuming translational invariance for the ground state as in Ref. (Q), 
we can write the master field as 



Upq = I J px+iqy U ah '{PA) 

Jn,n a 3 



'p,q a/3 

h*(p,«) = <ar)»*(p + 4 )-5 5 ^ rT — f , , ) 

o a/3 P 2 - m„ - P-r a+ (p) + 2e V 2p + - T a+ (p) -\J2m a ) t 

(2.1 
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Using the previous eq. ( |2.17| ) the eq. ( |2.16| ) can be easily rewritten in the 
momentum space as 

Tl(p) = 2g 2 J2(-i) n nf n i:(p,p) (2.18) 

n 

where 



JPL.-Pn-! (p ~ Pi)- (p 1 -p 2 )-" (Pn-1 

(2.19) 

where p_ is defined in eq. ( |2]7] ). 

Using the principal value regularization for t- as in Ref. ([[[1]]) and sym- 
metric p + integration done by 't Hooft (0, we get the result 



/oo 
dk x . . . dk n -i 
-oo 



V V V V 

r s 9n{kx)- —sgn(k 2 ) . . . sgn^kn^)- 

p - ki ki - k 2 k n - 2 - k n -.i k n -i - p' 

(2.20) 

Notice the useful property: I n (ap,ap') = ^I n (p,p'), that will turn out to be 
useful while computing the full qq amplitude in the next section. 



3 The full quark- ant i quark amplitude 

In order to compute the full qq amplitude we add a source J coupled to the 
composite field U to the action. The partition function in presence of an 
external source J is given by: 

Z[J) = Zj = J [dU][dM]e iNSj (3.1) 

where 

Sj = S eff - i Tr(JU) (3.2) 

Differentiating wrt J the logarithm of the partition function, we find the 
connected Green functions (from now on all the capital roman indices will 
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be used to indicate collection of indices (xaa)) 

1 8 log Zj 

Gab =< Uab >c= "tt 

1 N SJba 



j=o 



TT TT 1 S h & Z J 
GAB;CD —< UABUCD >c— TT^Tl FT 

2 N 2 5J B aSJdc 

To the leading order in N we have obviously that 

\og{zf)=iNSj[U] 

< U AB > { c ] = Uab 
o 

< UabUqd >c°^= JjUaB;CD 

where U is the saddle point value in presence of the source, i.e. 
j 

Uab = Uab + Uab-.pqJqp + ■ ■ ■ 
j o 1 

and it is given by the solution of the equation generalizing ( |2.14| ) in the 
presence of J 

2 

D + iU~ l + f - £ n/ n / (P( 2 :)C/) n - 1 P(a:) + J = (3.9) 

This equation cannot be solved exactly for an arbitrary J, but it can be 
solved perturbatively in J. Inserting ( [TH] ) in the previous equation and 
using 

(3.10) 



(3.3) 
(3.4) 

(3.5) 
(3.6) 

(3.7) 
(3.8) 



we find 



Uab ~ \U~ X -U^U-sqU^Jqp) 

J \0 1 ^ J AB 

2 

D +i zy^ + f-X: n/„ / (P(x)w) n ~ 1 P(x) = 

47T ~ Jx 



(3.11) 



and 



- iU~ l U; PQ U 
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-i 



iX.Qp + 

m=n—2 



& lib ^ p 

E /. 

^" n=2 m=0 17 x 



{P{x)U) m P{x) 



U-PQ 



[P(x)U) 



n—m—2 



P(x) 



■ 
(3.12) 
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where Xab,pq = SaqSbp The first of these equation has as solution the 
expression in eq. ( [2.17| ), while the other can be formally solved as follows^ 



m=n— 2 



u 

1 



\ ^" n=2 m=0 Jx 



n=2 

If we define 



(P(x)U) 



m+l 



n—m—l 



u®u l x 

o o 



iU®U l XQ 





m=n—2 
^" n=2 m=0 



(P(x)W) 



m+l 



:p(x)w) 



n— m— 1 



(3.13) 



(3.14) 



the previous equation for U can be rewritten in two different revealing forms 



as 



u 

1 



-{u ® - (w ® w*x) w 

1 



iU®U*)X + i{U®U t X) 





S( ^ + iu®u t xg)- 1 





W®W'-^.15) 





The first recursive form reveals clearly that the role of the usual gluonic 
interaction is now played by Q. Moreover the second form shows that 



r = G{ ^ + iXu t ®UG) 



-i 



(3.16) 



is the full quark-antiquark truncated amplitude and it is the generalization 
of the full qq amplitude of Callan, Curtis and Gross Tqcg (0)) as it is 
transparent from the following recursion relation, that it satisfies 



t = g-igx{u t ®u)T 





(3.17) 



Notice however that the recursion relation for T has a minus sign in front 
of the second term on rhs of difference wrt the CCG one: this has no conse- 
quences and it is due to the fact that we are dealing with an effective theory 



2 We use the convention that all the matrix operations written in capital letter refer to 

the mesonic indices, i.e. couples (AB) . . ., while the others to the quark ones, i.e. A, B 

Explicitly (A t )ab,cd = Acd,ab is the transposition wrt to the mesonic indices while 
(X )ab = Xba is wrt to the quark ones. 

We use also (X <E> Y)ab,cd = X ac Y B d and {A, PQ ® Y)ab,cd = A AC ;pqY B d and 
(A-pqX)ab = Aac-,pqXcb- 
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where the gluonic degrees of freedom are integrated away. This implies the 
appearance of quartic fermionic vertices and hence the substitution of every 
loop quark-gluon-quark-gluon in the original theory by a fermionic loop due 
to the elimination of the gluons. Since every fermionic loop carries a minus 
sign, this explains the origin of the difference. 

Now we go to the momentum space and we write the "gluonic" interaction 
Q and the full qq amplitude explicitly as 

<3ab-,cd = Q^s\p, q, q', P'} = (2n) 2 5\p + p' + q + q') 5 bc 5 da ( 7 _) />y (7_) fa G" 6 (p, p';r = p + q) 

(3.18) 

m=n— 1 

G ab [p,q,q' lP ']=27rg 2 £H) n n/„ £ I b m (q, -q')I^ m _ 1 (p' , -p) (3.19) 

n=2 m=l 

and 

Tab;CD = 7$%\p, q, q',p'] = (2n) 2 5 2 (p + p' + q + q') 5 bc 5 da M^h-) Sa T ab (p,p'; r) 

(3.20) 

where all the momenta are incoming. The recursion relation ( |3.17| ) reads 
now as 

T ab (p,p';r) = G ab (p,p';r) 

- 2i [ G ab (p,-m;r)p a _(-m)pt(r - m)T ab (m,p; r\3. 21) 



Defining as usual the x and x' variables as p- = xr_ , p'_ = — x'r_ (the minus 
is due to the fact that our momenta are all incoming and we want to make 
contact with the usual conventions) and 

<t>{x V x'-r) = \jdp + p a _(-p)p b (r-p)T ab (p,p';r) (3.22) 
the previous equation ( |3.21| ) yields T as a functional of (f) 

T ab (p,p';r) = ^G ab (x,-x';l)-^- [ dy (^(x^l)^,-^) (3.23) 
Defining the eigenfunctions $f (a?) through the equation 



M? ab)n <f>t(x,x';r) 



m a + xT + (x) m 2 + (1 — x)Y + {l — x) 
x 1 — x 



^{x.x'-r) 



1 dy -G ah (x, —y; l)^(y, x'\ r) (3.24) 

7T 
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and supposing that they form a complete system in j 6 [0, 1], we can write 
the full amplitude T as 

T ab (p.,p'_;r) = ^G ab (x,-x';l) 
En ^rzwr- Jo dy G ab (x, -y; r) j* dy' G ab (y', -x'; l)4>t*(y'; r) 



(3 



Specializing the previous equation to the normal case, we notice that Tqcg = 
2i T; the factor i can be traced back to the definition ( |3.16| ) and the factor 
2 is due to the definition of U ( |2]7] ) (see also ( |3/7| )). 



4 The general case. 

The most general expression for the generalized QCD2 is given by the action 

5 = Jd 2 x\ ^tr c (EF) - ^ E iW II M^) + m l 11)^ 

^ {n.i} ne{rii} 

(4.1) 

where {rii} is a whatever sequence of natural number satisfying the condition 
Hi < rij + i. Repeating the previous steps we find the effective action 

^S eff = Tr(DU+^-Y l f{m}J[ II (P(x)Ur]+MU-i\ogM 

\ {ni} X ne{rii} J 

(4.2) 

and the saddle point condition 
Sft+DEJW E jt P (x)U) n II Tr(P(x)U) m = -i 1 (4-3) 

{"J nG{n;} m £{ n .}\{ n } 

Using the translational invariance of the saddle point, i.e. of the ground 
state, it turns out that 

T n = Tr{P{x)U) n = ^(-^\ lY,IZ(p,P)p a -(p) (4-4) 

is independent of x, this means that the whole effect of introducing the most 
general potential amounts in generating the effective coupling constants /„ 
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given by 

In = Y E f{n i} d n ({ni}) II T ™ (4-5) 

n {ni} such that ne{m} me{nj 

where d n ({rii\) is the degeneration of n in the sequence {rij}. As it can be 
shown easily all the formula of the previous section are still valid with the 
substitution f n — > f n (|4.5| ). 

We notice however that the theory is independent of f% because J x P(x) = 
as a consequence of J x P(x) oc a+7- and of the absence of a constant vector 
a + in the theory. 

Moreover we have 

T 2n+1 = (4.6) 

and this implies that it is impossible to induce terms with odd powers starting 
from terms with even total powers. The demonstration of (|4.6| ) goes as 
follows: 

Tn + 1 OC / 1 ffi?(pih 1 -^T ■ ■ ■ (£ b (Pn-l) , 1 v P b -( P ) 

Jp,pi...p n [P ~ Pi)- [Pi ~ P2)- [Pn ~ P)- 

(4.7) 

then we rename the integration variables as p <-> p n ,Pi —> p n -i,P2 ~^ Pn-2 ■ ■ ■ 
and we use p~(— p) = — (that follows from fl2.18| , |2.17j )), and we find 

T n = (— ) n T n , that proves what asserted, independently of the regularization 
scheme. 



5 Conclusions. 

In this work we have solved the most general case of minimally coupled 
generalized QCD 2 (HTT||) using the bilocal method ((jTJ, ||). This is achieved 



in a straightforward and transparent way, thus showing the effectiveness of 
the method, when it is applied to vector-like models coupled to matrix gluon 
fields, which in two dimensions are not dynamical (see also [|13|| ): the main 



point that guarantees the successful application of the large N techniques 



is the possibility of defining colourless fields (see eq. (|2.7|) ) that allows 
one to extract the iV dependence both in the action and in the measure of 
integration (the mesonic bilocal U(x, y) is a global colour singlet, as it is the 
auxiliary bilocal field M(x,y)). 
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